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Abstract 

For  over  10  years  primary  timing  laboratories  have  used  signals  from  the  Global  Positioning 
System  (GPS)  to  compare  atomic  timescales.  The  method  used,  the  common-view  of  GPS  satellites, 
does  not  use  the  GPS  system  to  its  full  potential.  Recently,  there  has  been  much  interest  in  using 
geodetic-quality  GPS  receivers  for  time  transfer.  The  result  has  been  a  substantial  improvement  in 
the  precision  of  the  resulting  time  and  frequency  transfers. 

In  this  paper,  a  detailed  account  is  presented  of  common-clock  comparisons  made  at  NPL  between 
a  recently  purchased  Ashtech  Z12-T  receiver  and  an  older  Allen  Osborne  Associates  (AOA)  TTR-4P 
receiver.  Data  collected  from  these  geodetic-quality  GPS  receivers  were  processed  using  analytical 
software  developed  at  NPL.  Direct  comparisons  are  made  between  the  two  receivers  using  both  PI 
and  P2  coded  signals,  and  LI  and  L2  phase  measurements.  Measured  ionospheric  delays  were 
obtained  from  both  receivers  and  compared.  Multipath  is  one  of  the  major  causes  of  errors  in  the 
surveying  applications  of  geodetic  GPS  receivers.  Code-phase  differences  were  used  to  estimate  the 
magnitude  of  code  multipath  present  in  both  receivers.  The  principal  sources  of  errorpresent  in  the 
common  clock  measurements  are  discussed  and  possible  improvements  are  considered.  The  future 
direction  of  geodetic  time  transfer  work  at  NPL  is  also  outlined. 


1  INTRODUCTION 

The  application  of  geodetic  GPS  techniques  to  frequency  and  time  transfer  has 
produced  some  very  impressive  results  [l],[2j.  These  results  have  increased  the 
interest  in  using  geodetic  GPS  time  transfer  for  international  comparisons  between 
primary  timing  laboratories.  One  recent  development  has  been  the  formation  of  a 
joint  B1PM/IGS  pilot  project  to  coordinate  studies  in  the  field  [3]. 

The  quality  of  geodetic  GPS  frequency  and  time  transfers  will,  however;  be  limited  by 
the  delay  stability  of  both  the  GPS  receiver’s  instrumentation  and  that  of  the 
associated  measurement  systems  [4].  Other  local  environmental  considerations,  for 
example  the  susceptibility  of  the  receiver’s  antenna  to  multipath,  will  also  limit  the 
receiver’s  performance. 
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The  post-processing  of  geodetic  GPS  data  is  more  complicated  than  the  processing  of 
data  collected  by  other  high  precision  time  transfer  methods.  Only  a  few  primary 
timing  laboratories  have  the  capability  to  postprocess  their  own  geodetic  GPS  data, 
and  those  who  do  usually  use  commercial  software. 

In  this  paper  the  results  of  a  very-short-baseline  common-clock  experiment  between 
NPL’s  Ashtech  Z12-T  and  AOA  TTR-4P  geodetic -quality  GPS  receivers  are 
presented.  NPL’s  development  of  analysis  software  from  first  principles  is  described. 
The  relative  performance  of  NPL’s  GPS  receivers  is  evaluated  by  using  the  analysis 
software  to  compare  both  the  code  and  phase  data  available  from  both  receivers.  By 
combiningthesedata  with  the  broadcast  GPS  almanac  data,  a  detailed  evaluation  was 
possible,  including  an  examination  of  the  physical  processes  limiting  the  receivers’ 
performance. 

2  GEODETIC  GPS  HARDWARE  AT  NPL 

NPL  possesses  two  geodetic  quality  GPS  receivers,  an  Ashtech  Z12-T  and  an  AOA 
TTR-4P.  The  hardware  configuration  is  shown  in  Figure  1.  The  Ashtech  Z12-T  is  a 
modified  version  of  the  Ashtech  Z12  receiver,  where  the  internal  oscillator  is 
replaced  by  a  20  MHz  external  reference  signal.  The  receiver  provides  1  PPS  input 
and  output  signals,  which  lock  the  receiver’s  internal  clock.  The  TTR-4P  is  an  older 
receiver  widely  used  in  primary  timing  laboratories.  The  two  receivers’  antennas  are 
located  approximately  five  meters  apart.  Standard  frequency  and  1  PPS  reference 
signals  for  the  two  receivers  originate  from  common  frequency  distribution  and  1 
PPS  generation  units.  All  the  receivers’  input  signals  are  referenced  to  NPL’s  Sigma 
Tau  Active  Hydrogen  Maser. 

3  DATA  SETS  EXTRACTED  FROM  THE  RECEIVERS 

Several  data  sets  were  extracted  from  the  receivers.  RINEX  data  were  extracted  from 
the  Ashtech  Z12-T  receiver.  Although  in  principle  available,  similar  RINEX  data 
could  not  be  obtained  from  the  TTR-4P  receiver;  however,  code  and  phase 
measurements  were  obtained  from  the  receiver  using  its  “block  015"  data.  Both  data 
sets  were  collected  with  a  measurement  epoch  of  30  seconds  synchronized  to  GPS 
time.  The  GPS  almanac  data  were  also  extracted  from  the  Ashtech  Z12-T  receiver, 
which  was  used  in  NPL’s  postprocessing  software. 

4  POSTPROCESSING  SOFTWARE  DEVELOPED  AT 
NPL 

The  baseline  between  NPL’s  two  geodetic  GPS  receivers  is  5  m.  With  this  limited 
baseline  extensive  data  postprocessing  is  still  required  to  compute  comparisons 
between  the  two  receivers  with  a  measurement  precision  of  10  picoseconds. 
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Several  distinct  operations  are  performed  by  the  analysis  software: 

1)  Each  day’s  data  are  processed  separately.  The  software’s  data  input  included  code 
and  phase  variables  from  each  receiver,  along  with  the  GPS  almanac  data  and 
approximate  receiver  coordinates.  The  data  are  processed  in  30  s  blocks,  which 
corresponds  to  the  measurement  epoch. 

2)  The  GPS  almanac  is  decoded  to  obtain  satellite  coordinates,  at  each  measurement 
epoch,  which  were  then  used  to  calculate  accurate  satellite  azimuth  and  elevation 
values.  Satellite  azimuth  and  elevation  angles  are  generated  in  the  Ashtech  Z12-T 
and  TTR-4P  output  data  sets,  but  only  with  a  precision  of  1  °,  which  is  not  adequate 
for  the  comparison  of  the  receivers’  clocks  with  10-picosecond  precision. 

3)  A  first  calculation  is  made  of  the  relative  coordinates  between  the  phase  centers  of 
the  two  antennas,  using  only  the  code  variables  PI  and  P2,  and  the  satellite  azimuth 
and  elevation  values  calculated  in  step  2.  A  separate  estimate  is  made  at  each 
measurement  epoch  by  calculating  the  pseudo-range  differences  (P1<ttr-4P)-P1  (Ashtech 
zi2-T))  and  (P2(TTf.4p)-P2(Ashtech  zi2-T))  using  each  satellite.  The  antennas’  relative 
coordinates  are  then  optimized  to  minimize  the  scatter  between  pseudo-range 
differences  calculated  from  each  satellite  in  view.  This  calculation  is  repeated  for 
each  measurement  block,  and  finally  a  mean  value  for  the  relative  coordinates  is 
determined. 

4)  Phase  dataarenow  introduced  into  the  analysis.  Step  (3)  is  repeated  using  phase 
rather  than  code  measurements.  Phase  ambiguities  are  corrected  by  adding  integer 
cycles  to  the  LI  and  L2  measurements.  With  these  very  short  baselines,  resolving  the 
phase  ambiguities  is  relatively  straightforward.  Using  a  single  day’s  data,  the 
coordinate  differences  between  the  two  antennas  could  be  determined  with 
sub-centimeter  precision.  Separate  relative  coordinates  are  calculated  for  the  LI  and 
L2  frequencies  to  reflect  the  slightly  different  phase  centers  in  the  antennas. 

5)  The  software  is  then  used  to  generate  a  wide  range  of  output  parameters  using  the 
optimized  relative  coordinates  in  the  processing.  These  include:  measured  code  and 
phase  differences  between  the  two  receivers,  ionospheric  delay  measurements,  and 
receiver  code/phase  differences  used  in  the  study  of  code  multipath. 

5  DIRECT  COMPARISONS  BETWEEN  THE  TWO 
RECEIVERS 

PI  code  differences,  calculated  between  the  two  receivers  are  shown  in  Figures  2  and 
3.  P2  code  differences  were  also  calculated  but  are  not  shown  here.  Figure  2  shows 
individual  satellite  tracks,  while  the  mean  value  calculated  from  all  satellites 
simultaneously  tracked  by  both  receivers  is  shown  in  Figure  3.  The  PI  and  P2  code 
differences  showed  no  obvious  correlations.  There  is  a  scatter  of  several  nanoseconds 
within  individual  tracks:  this  decreases  to  (1-2)  nanoseconds  (lo)  for  the  mean  value 
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averaged  over  all  satellites.  The  scatter  in  the  data  appeared  asymmetric,  with  notable 
outliers.  The  majority  of  these  were  attributed  to  multipath  occurring  within  the  AOA 
TTR-4P  receiver,  often  when  the  satellite  is  at  low  elevation  angles.  The  data  scatter 
reduces  with  satellite  elevation  angle,  due  in  part  to  the  increased  gain  of  the  antenna 
at  these  angles.  There  was  no  obvious  correlation  between  individual  satellite  tracks. 
The  scatter  on  the  P2  data  is  approximately  twice  that  on  the  PI  data.  This  is 
consistent  with  the  lower  levels  of  received  P2  signals. 

The  phase  differences  between  the  two  receivers  are  shown  in  Figures  4  and  5.  The 
scatter  on  the  phase  data  is  noticeably  lower  than  that  from  the  code  data,  typically 
1 0-picosecond  scatter  between  successive  measurements  on  the  same  satellite  track . 
The  long-term  variations  show  no  obvious  correlation  with  long-term  variations  in 
the  code  measurements.  Individual  satellite  tracks  correlated  extremely  well  (Figure 
4), with  a  typical  20-picosecond  RMS  scatter  (la)  between  all  the  tracks  measured  in 
the  same  epoch.  This  shows  that  the  majority  of  the  phase  instabilities  are  correlated 
between  measurements  from  all  of  the  satellites.  The  phase  differences  observed  at 
LI  and  L2  frequencies  also  agreed  well  (Figures  5), although  there  were  variations  of 
up  to  60  picoseconds  in  the  L1-L2  phase  differences  (Figure  6).  The  short-term 
instabilities  in  the  phase  differences  correlate  with  temperature  cycling  in  the  room 
containing  the  receivers;  this  cycling  has  a  period  of  approximately  30  minutes  (See 
Figures  4,5,  and  6).  The  Ashtech  Z12-T  receiver  temperature  is  included  in  Figure  5. 
There  was  no  obvious  correlation  with  outdoor  temperature.  More  detailed  studies 
are  planned  to  isolate  the  causes  of  the  phase  instabilities. 

6  IONOSPHERIC  DELAY  MEASUREMENTS 

Over  a  longer  baseline  a  geodetic  GPS  time  transfer  will  be  limited  by  the 
effectiveness  of  the  ionospheric  delay  correction.  Three  different  methods  of 
measuring  the  LI  ionospheric  correction  are  shown  here:  (i)  using  code  a(Pl-P2) 
differences  (Figure  7);  (ii)  using  phase  a(L2-Ll)  differences  (Figure  8);  and  (iii) 
using  single  frequency  code-phase  0.5(P1-L1)  differences  (Figure  9)  calculated  at 
the  LI  frequency,  (a  =  F22/(F22-Fi2)  where  Fi  and  F2  are  the  transmission  frequencies 
of  LI  and  L2  respectively).  Absolute  values  of  the  ionospheric  delay  may  be 
obtained  from  the  code  measurements;  however,  only  relative  ionospheric  delay 
changes  may  be  obtained  using  the  phase  measurements,  and  the  code-phase 
differences.  Good  agreement  is  obtained  between  all  three  methods.  The  phase 
measurements  clearly  demonstrated  far  less  noise  than  on  the  code  measurements; 
this  was  in  part  due  to  the  encryption  of  the  GPS  P-code,  and  in  part  due  to  phase 
measurements  being  intrinsically  less  noisy.  The  noise  on  the  code  measurements 
obtained  from  the  AOA  TTR-4P  receiver  was  approximately  half  the  level  of  noise 
obtained  from  the  Ashtech  Z12-T  receiver.  There  was  no  obvious  bias  between  code 
measurements  obtained  from  different  receiver  channels  or  from  different  GPS 
satellites;  however,  there  was  a  bias  of  30  nanoseconds  between  the  coded 
ionospheric  measurements  obtained  from  the  Ashtech  Z12-T  and  AOA  TTR-4P 
receivers.  The  single  frequency  code-phase  differences  are  shown  to  be  effective  in 
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determining  ionospheric  delay  changes. 

The  differences  between  ionospheric  phase  measurement  obtained  from  the  Ashtech 
and  AO  A  TTR-4P  receivers  are  shown  in  Figure  10.  Variations  in  these  plots  of  up  to 
1 00  picoseconds  were  observed  which  correlate  with  satellite  elevation.  This  may  be 
due  to  the  different  relative  phase  patterns  of  the  two  antennas.  A  simple  alterative 
explanation  is  that  the  height  component  of  the  antenna  coordinate  differences  may 
be  slightly  in  error. 

7  MEASUREMENT  OF  CODE  MULTIPATH 

A  major  limitation  to  the  use  of  geodetic  GPS  receivers  for  positioning  applications 
is  multipath.  Phase  multipath  is  typically  a  thousand  times  smaller  than  code 
multipath.  The  quality  of  the  GPS  code  measurements  is  vital  for  the  analysis  of 
geodetic  GPS  data,  which  may  be  severely  degraded  by  code  multipath.  Measuring 
code-phase  differences  from  individual  satellite  tracks  is  an  excellent  way  to  identify 
code  multipath.  The  repeating  GPS  constellation  each  sidereal  day  enables  multipath 
to  be  observed  through  the  repeating  delay  changes  observed  in  the  code  phase 
differences.  Examples  are  shown  in  Figure  1 1  for  data  obtained  from  NPL’s  Ashtech 
Z12-T  receiver.  The  repeating  code  phase  difference  patterns  were  observed  at  both 
LI  and  L2  frequencies.  The  patterns  observed  were  noticeably  different  at  the  two 
frequencies,  because  multipath  effects  are  frequency-dependent  (Figures  1 1  and  1 2). 
The  origin  of  much  of  this  multipath  may  well  be  a  copper  dome  on  NPL’s  laboratory 
roof.  The  Ashtech  Z12-T  antenna  has  recently  been  moved;  the  new  site  will  be 
examined  for  residual  code  multipath. 

8  FUTURE  WORK 

NPL  has  an  active  program  to  develop  its  geodetic  GPS  installation.  NPL  intends 
to  purchase  another  geodetic-quality  GPS  'receiver,  enabling  “  three-cornered-hat” 
comparisons  in  future.  Also  NPL  has  a  3S  Navigation  combined  GLONASS/GPS 
receiver  from  which  the  LI  CA  code  GPS  channels  may  be  used  for  direct 
comparisons  with  geodetic  GPS  receivers.  Instrumentation  delay  instabilities  may  be 
reduced  through  temperature  control  of  the  GPS  receivers,  antennas,  and  connecting 
cables.  Work  is  underway  to  improve  the  frequency  distribution  system  from  NPL’s 
Active  Hydrogen  Maser. 

NPL  is  in  the  process  of  upgrading  its  geodetic  GPS  installation  to  meet  the 
requirements  of  an  IGS  station.  NPL  is  also  playing  an  active  part  in  comparisons 
between  geodetic  GPS  and  other  high  precision  time  transfer  methods,  including 
TWSTFT  and  GLONASS  time  transfer. 
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9  CONCLUSIONS 


This  paper  has  demonstrated  that  short- baseline  common-clock  experiment  using 
geodetic  GPS  receivers  may  be  analyzed  from  first  principles  using  relatively 
straightforward  software.  This  type  of  analysis  may  be  usefully  performed  by  primary 
timing  laboratories  to  obtain  an  insight  into  the  performance  of  their  geodetic  GPS 
timing  receivers. 

A  major  source  of  errors  is  due  to  delay  instabilities  within  the  receiver 
instrumentation.  The  analysis  described  in  this  paper  has  identified  the  magnitude  of 
both  the  code  and  phase  instabilities  occurring  between  NPL’s  geodetic  GPS  timing 
receivers.  Other  physical  processes  limiting  geodetic  GPS  time  transfer,  including 
ionospheric  delay  determination,  code  multipath,  and  antenna  phase  dispersion,  have 
been  examined.  Future  work  is  required  to  optimize  the  receiver’s  performance. 
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MJD-50913 


Figure  2:  PI  code  comparisons  (AO A  TTR-4P  -  Ashtech  Z12-T),  individual  satellite 
tracks. 
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Figure  7:  Code  ionospheric  delay  measurements,  individual  satellite  tracks,  Ashtech 
Z12-T  receiver. 


Figure  8:  Phase  ionospheric  delay  measurements,  individual  satellite  tracks.  Ashtech 
Z12-T  receiver. 
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Questions  and  Answers 


MIHRAN  MIRANIAN  (USNO):  On  the  TTR4-P.  were  you  using  any  kind  of  an  elevation  mask  for  the 
ionosphere?  Because  it  is  notoriously  bad  for  ionospheric  measurements. 

JOHN  DAVIS  (NPL):  No,  we  were  not.  We  were  using  everything  that  we  could  collect  from  the  TTR4-P 
because  the  Ashtech  had  more  channels  (12);  the  TTR4-P  had  eight.  Also,  we  discovered  that  the  TTR4- 
P  found  it  slightly  harder  to  lock  onto  satellites.  So  in  reality,  everything  was  used.  Where  there  was  a 
TTR4-P  track,  there  was  always  an  Ashtech.  We  did  not  use  a  mask.  In  those  P-1  measurements, 
everything  has  been  included,  wherein  after  the  average  including  the  low-elevation  satellites  was 
measured,  we  could  put  a  mask  on. 

MIHRAN  MIRANIAN;  Yes,  we  use  them  at  the  Observatory,  and  we  found  that  anything  below  about  25 
degrees  is  realty  useless. 
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